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Entanglement catalysis allows one to convert certain entangled states into others by the temporary 
involvement of another entangled state (so-called catalyst), where after the conversion the catalyst is 
returned to the same state. For bipartite pure entangled states that can be transformed in this way 
with unit probability, the respective Schmidt coefficients are said to satisfy the trumping relation, 
a mathematical relation which is an extension of the majorization relation. This article provides all 
necessary and sufficient conditions for the trumping relation in terms of the Schmidt coefficients. 
The coefficients should satisfy strict inequalities for the entropy of entanglement and for power 
means excluding the special power f . 
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I. INTRODUCTION 

An important problem in quantum information the- 
ory is to understand the conditions for transforming a 
given entangled state into another desired state by us- 
ing only local quantum operations assisted with classi- 
cal communication (LOCC). Significant development has 
been achieved for the case of pure bipartite states. In the 
asymptotic limit, where an infinite number of copies of 
entangled pairs are to be transformed into each other, it 
is found that the conversion is possible with the proba- 
bility of success approaching unity as long as the entropy 
of entanglement does not increase. [l[ 

However, away from the asymptotic limit, where a sin- 
gle copy of a given state is to be transformed into an- 
other given state, such a simple conversion criterion can- 
not be found and investigations have unearthed a deep 
connection of the problem to the mathematical theory 
of majorization. 2] For two sequences of n real numbers 
X = {xi,X2, ■ ■ ■ ,Xn) and y = (2/1,2/2, ■ • • ,?/«), we say that 
X is majorized by y (written x ^ y) when 

xl+xl + --- + xl>yl+yl + --- + yl (1) 

for m — 1,2, . . . ,n — 1 and the sequences have the same 
sum — Til Hi)- Here, x^ represents the sequence x 

when its elements are arranged in non-decreasing order 
{x\ < x\ < ■ ■ ■ < x"}^) and similarly for y^ . Nielsen has 
shown that for two given states having the Schmidt forms 

n 

IV") = ^^Jxl\^A®^B) , (2) 

i=l 
n 

\<l>) = , (3) 

where x and y are the respective Schmidt coefficients 
= = 1)1 the state IV') can be converted 

into by LOCC with unit probability of success if and 
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only if X ^ y-W\ Subsequently, Vidal extended this result 
to probabilistic transformations where the probability of 
success is related to the violation of the majorization in- 
equalities in (dJ.Q 

Soon afterwards, Jonathan and Plenio have demon- 
strated an interesting effect that is termed as catalysis 
or entanglement assisted local transformation, [sj There 
are some cases where j?/') cannot be converted into \(j)) 
with certainty (conversion is possible only with a proba- 
bility less than 1), but with the involvement of another 
entangled pair (a catalyst), the conversion is made pos- 
sible. In other words, if |x) = TIi=i \fci\^A ® ^b) is the 
state of the catalyst, then |V') ^ lx) can be converted into 
10) ® Ix) with complete success. In such a transforma- 
tion, the entanglement of the catalyst is not consumed, 
although it takes part in the transformation. Catalysis 
is also useful in almost all conversion processes where it 
improves upon the conversion probability. 

Expressing in terms of the Schmidt coefficients and 
considering only the cases where catalysis helps achieve 
unit probability of success, we basically have situations 
where x is not majorized by y, but there is a sequence c 
such that x®c\% majorized by y (8) c. Following Nielsen, 
for two sequences of non-negative numbers x and y with 
n-elements, we will say that x is trumped by y (written 
X ~<T y) if there exists another sequence c such that x ® 
c < y^c. It is easy to see that in such cases the catalyst 
sequence c can be chosen from strictly positive numbers. 

A lot of research has been directed to understand the 
catalytic transformations @ and to analyze the mathe- 
matical structure of the trumping relation. 0, H, H, 
One of the open problems is to find a way to decide if 
two given sequences x and y satisfies the trumping rela- 
tion. The purpose of this article is to give all necessary 
and sufficient conditions for this relation. This problem 
has been partially solved by Aubrun and NechitapTl.[T3|. 
who work with stochastic tools to describe the closure of a 
set constructed with the trumping relation. The methods 
used in this article, however, are quite similar to those of 
a recent study that provided an expression for the cat- 
alytic conversion probability. [l^ But, as the trumping 
relation necessarily implies that the two sequences have 



2 



the same sum, the mathematical details of the proofs 
given below are more complicated than the ones in that 
article. 

It appears that the necessary and sufficient conditions 
for the trumping relation can be expressed in terms of 
the power means, which are defined as 



and the entropy of entanglement, 

n 

a{x) — — Xi \n c 



(4) 



(5) 



The main theorem that we would like to prove is the 
following. 

Theorem 1. For two n-element sequences of non- 
negative numbers x and y such that x has non-zero el- 



ements and the sequences are distinct (i.e. 



the relation a; ?/ is equivalent to the following three 
inequalities 



A^{x) < A^{y) 
(t{x) > a{y) 



y i^e (-00,1) 
V I' e (l,+oo) 



(6) 
(7) 
(8) 



where all inequalities are strict. 

Note that by the continuity of the power mean func- 
tion against v, the requirement that the sequences x and 
y have the same sum is included in the conditions ^ and 
([7]). Moreover, the limits of these inequalities at = — cxd 
and v = +00 imply that the minimum and maximum ele- 
ments of the sequences satisfy the respective inequalities 
x\ > yl and a;Jj < yj,, but these do not have to be strict. 
Note also that all of these three conditions can be equiv- 
alently expressed as the strict positivity of the function 



1 



1 



■In 



A.jy) 
A^{x) 



(9) 



for all finite values of i^, where the value Ri corresponds 
to the difference of the entropies. 

Notice that some of the results in Ref. [3 about 
the connection between the multiple-copy entanglement 
transformation, a related phenomenon discovered by 
Badyopadhyay et aLp^. and the trumping relation can 
be easily understood in view of the conditions ([5][5]). 
For any sequence x and any integer fc > 1, we have 
= A^ix)'' and cr(a;®'=) = ka{x). As a result, 
if k copies of a state with coefficients x can be trans- 
formed into k copies of another state with coefficients y, 
either with or without catalysis, then x must be trumped 
by y. 

The article is organized as follows: In section |lll a 
few relations related to majorization are given and a key 
lemma is proved. Then, in section IIIIl the theorem is 
proved. Finally, section HVl contains the conclusions. 



II. A FEW RELATIONS AND A KEY LEMMA 

The following facts about the majorization and the 
trumping relation will be used occasionally. 

(1) If X ^ y, then for any convex function /, we have 

n n 
1=1 i=l 

Moreover, if x^ ^ y^ and / is strictly convex, then 
the inequality above is strict. 

(2) For any sequence a;, we define the characteristic 
function Hx{t) = X]r=i(^ ~ ^i)^ where (a)+ ~ 
max(a, 0) denotes the positive-part function. For 
non-negative sequences x and y which have the 
same sum = S^i); the relation x -< y can 
be equivalently stated as 



Hx{t)<Hy{t) Vi>0 



(11) 



(3) For the cross-product of two sequences we have 
Hxi^c — '}2i'^f-Hx{t/ce). 

(4) Relation -< and -<t are partial orders on all n- 
element sequences (up to equivalence under rear- 
rangement). Moreover, x ^ y implies x -<t y- 

(5) Let z be any sequence of strictly positive numbers 
and let x(B z denote the sequence obtained by con- 
catenating the elements of z to those of x. Then 
X -<T y X® z y®z. Moreover, x and y satis- 
fies the inequalities ([6][11) iff a; z and y®z satisfies 
the same inequalities. As a result, nothing will be 
lost from generality if a proof of the Theorem 1 is 
given for sequences that have no common elements. 

The proof of the sufficiency of the conditions is 
based on the following key lemma. 

Lemma: If a polynomial 7(5) has no positive roots 
and 7(0) > then 

(a) it can be expressed as 7(5) = b{s)/a{s) where a{s) 
and b{s) are polynomials with non-negative coeffi- 
cients. 

(b) Moreover, a(s) can be chosen as a polynomial with 
integer coefficients. 

Proof: The lemma can easily be generalized to poly- 
nomials which have a root at s = 0, but for our purposes 
the above form is sufficient. For part (a), we first provide 
the proof for a second degree polynomial with complex 
roots, e.g., 7(5) = 1 — 2^3 + As^ where A > £_^. For 
^ < 0, there is nothing to be proved as 7 has already 
non-negative coefficients. For ^ > 0, 7(5) can be written 
as the ratio 6(s)/a(s) where 



2N-1 



- E(l + A.2)'=(2e.; 



fc=0 



is) = (H-A.^) 



2\2N 



(13) 



2N 
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and iV is a sufficiently large integer so that the inequality 



1 /(2jv)!\^ e 



(14) 



is satisfied. It is always possible to find such an N since 
the left-hand side has a limit 1 when TV goes to infinity 
and the right-hand side is strictly less than 1. For such 
a choice of N, both a{s) and b{s) will have non-negative 
coefficients. 

For a general polynomial 7 which has no positive root, 
we first express it as a product of its irreducible factors 
as 



7(s) =A]1{1 + Qs) - 2C^s -I- i^f + 



(15) 



where A > 0, — 1/Ci are the real roots of 7 (therefore 
Q > 0) and {S^i±irii)~^ are the complex roots of 7. As the 
quadratic factors can be expressed as a ratio and the rest 
is simply a polynomial with non-negative coefficients, 7 
can be expressed as a ratio of two polynomials with non- 
negative coefficients. Note that, as 7 has no root at 0, 
both a{s) and 5(s) can be chosen as polynomials having 
a non-zero constant term. 

Before passing on to the proof of the statement (b), 
we first show that the polynomial b(s) can always be 
chosen such that all of its coefficients are strictly posi- 
tive. For this, consider a degree m solution for b{s), i.e., 
K^) = SfcLo ^ks'' where 60 > 0, b„i > and bk > for all 



1 < k < m. Let e(s) = 1 -I- s 



Then e{s)b{s) 



is a polynomial with degree 2m — 1 and all of its 2m coef- 
ficients are positive. Moreover, the polynomials e(s)6(s) 
and e(s)a(s) satisfy the conditions of part (a). Therefore, 
b{s) can be chosen to have non-zero coefficients. 

For the proof of part (b), suppose that b{s) is a 
degree m polynomial with positive coefficients and let 
(3 ~ mino<fc<m bk be the minimum of those. Let 



f3 



(16) 



where 7„ are the coefficients of the polynomial 7(s). De- 
fine a new polynomial a(s) such that it has the same de- 
gree as a(s) and its coefficients are chosen from rational 
numbers such that 



\ak - ak\ <e 0, 1, . 



.N 



(17) 



where and are the coefficients of a{s) and a{s) 
respectively. As the rational numbers are dense, this can 
always be done. If 0(5)7(5) = b{s), then the coefficients 
b{s) satisfy 



E 



'i\>~P (.18) 



Therefore, bk > bk — (3 > 0, i.e., b{s) has non-negative 
coefficients as desired. Multiplying a{s) by the common 
denominator of its coefficients gives a polynomial with 
integer coefficients. □ 



III. PROOF OF THEOREM 1 

Proof of the necessity of the conditions ([Bllll) for the 
trumping relation is trivial. Given that there is a catalyst 
c so that we have x(^c -< y(^c, we use the strict inequality 
of (flOl) for the following strictly convex functions: f{t) = 
f for > 1 and < 0, f{t) = for < t < 1, 
f{t) = -Int and f{t) = tint. All inequalities piSjl 
follow from these. 

The proof of the sufficiency of the conditions ([6][8]) is 
lengthy and requires us to separate it into three special 
cases. The key proof is for case A, where only the se- 
quences which can be expressed as integer powers of a 
common number is considered. Case B concentrates on 
non-zero sequences and uses the stability of the suffi- 
ciency conditions under small changes of sequences to 
reduce the problem to case A. Finally, case C deals with 
the situation where y has zero elements. 

Case A. y has strictly positive elements such that yi = 
Klli"-' and Xi — Kui^' for some integers ai and f3i and for 
some numbers K > and a; > 1. 

Proof: Without loss of generality, it is assumed that x 
and y have no common elements and they are arranged 
in non-decreasing order. The smallest of the exponents 
is ai which can be set equal to by a redefinition of 
K. Finally, both x and y can be divided by K which is 
equivalent to setting K = 1. As a result, it is not required 
that the sequences are normalized (i.e., they do not add 
up to 1). Since ai = 0, all other exponents satisfy ai > 
and Pi > 0. 

Let the polynomial r(s) be defined as 



kS 



(19) 



First, note that r(s) has simple roots at s = 1 and s — uj. 
This can be simply seen by evaluating its derivative at 
these points. 



r'(l) = ^(a,-/3,)<0 
1=1 

r'(..) . > 

In a; 



(20) 
(21) 



where the former strict inequality follows from ([6]) at = 
and the latter follows from ([5]). The fact that x and y 
are not normalized do not invalidate the latter inequality. 

Therefore, 7(3) = r(s)/((l — s)(l — s/ui)) is a polyno- 
mial. It can be seen that 7(0) > 0. Moreover, we will 
show that 7(5) has no positive root. For this purpose let 
s = uj'^ where v is any real number {v = and v — I can 
be excluded if desired). Then 



7K) 



iy- - 



(22) 



which can be seen to be strictly positive by virtue of ^ 
and ([7]) for all values of v. (For v — Q and v — 1, we have 
seen above that 7(5) has no root at 1 and lo). 
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By the lemma, there are two polynomials a(s) and b(s) 
with non-negative coefficients such that 0(5)7(5) — 6(5) 
and 0(5) has integral coefficients. The constant coeffi- 
cients a(0) and 6(0) will also be chosen to be non-zero. 
In terms of F, the relation can be expressed as 



a(5)r(5) = (1 - 5)(1 - s/uj)b{s) 



(23) 



Let a{s) has degree N. The catalyst sequence c will 
be chosen from the numbers w'^ which are repeated Uk 
times {k = 0, 1, ... , N). In that case, the characteristic 
function of c is 



N 



(24) 



fc=0 



We would like to show that the function 

A(t) = Hy^,{t) - H,^,{t) (25) 



1=1 

E 
E 



yiHc(t/yi) - XiHc(t/xi) 



is non- negative for alH > 0. First note that 



M+l 

E 

771=0 



(/m-/™-i)(t-u;")' 



(26) 
(27) 
(28) 



(29) 



where /(5) = (1 — s/uj)b{s), fm are coefficients of the 
polynomial /(s) and we have chosen /_i = for simplic- 
ity. Here M is the degree of / (M + 1 is the degree of 
a(5)F(5)). Since A(t) is a piecewise linear function, for 
showing its positivity, it is sufficient to look at its value 
at the turning points and at the and 00 limits. First 
note that A{t) = for t < 1 and A(t) is constant for 
t > As a result, we only need to check the values 



of A{t) ai t — Lu,L 



, cj*^+i. For any 1 < fc < Af + 1, 



fe-i 



A(u;'^-) = ^(/™-/„_i)(c.'=-c.") (30) 

fc-1 fe-1 
= {cJ-l)J2ifm-frn-l)J2''' (31) 



m=Q 

fc-1 p 



p—m 



(^-l)^a>P^(/™-/^_i) (32) 

p^O m— 

k-1 



(33) 



p=0 



Finally, as /(5) = (1 — 5/^)6(5), the coefficients of these 
polynomials satisfy 



fp 



Op-l 



(34) 



where 6_i = 0, which leads to 

A{uj'') = (w - l)6fc_iw'^-i > 



(35) 



This completes the proof of A{t) > for all t > 0. It also 
shows that x ® c < y ® c. Therefore, x -<t JZ-D 

Before passing on to the next case, we first state an- 
other theorem that shows the stability of the inequalities 
(OS]) against small variations in sequences x and y. Since 
only sequences with non-zero elements will be considered 
in the next case, the distance between two sequences will 
be measured by the deviation of the ratio of the corre- 
sponding elements from 1. For two sequences x and x 
which has no zero elements, the distance between them 
is defined as 



D(x] x) 



In- 



(36) 



The following theorem expresses the stability of the con- 
ditions (HlHl). 

Theorem 2. Let x and y be n-element sequences 
formed from positive numbers such that x[ > yj and 
xjj < y\. If X and y satisfy the inequalities ([6|), ^ and 

, then there is a positive number e such that whenever 
D{x\x) < e and D{y;y) < e, and Y,Xi = J^Vi = J2^i^ 
the sequences x and y satisfy the same strict inequalities. 

The proof of Theorem 2 is postponed to Appendix [X] 
This result will be used in the proof of the next case. 

Case B. y has strictly positive elements. 

The proof will be carried out by choosing two new 
sequences x and y which are sufficiently near to x and y 
such that Theorem 2 can be invoked, and it will be made 
sure that x and y satisfy the conditions considered in case 
A. Without loss of generality, it is assumed that x and y 
are normalized (^ Xi — yi — ^) and they are arranged 
in non-decreasing order. Let H = a{x) — a{y) > be the 
entropy difference of these sequences and let L = | In y j | . 
Note that the logarithm of all elements are bounded by 
L, i.e., |ln?/i| < L and llnx^j < L. Let eo be a positive 
number such that whenever D{x,x) < eo and D{y,y) < 
eo, the sequences x and y satisfy all the inequalities in 
(OS]). The positive number e is chosen such that 



^ . , eo 1 1 H \ 

e < mm — , — , — , 

2 ' 8n' n2' 96nL J 



(37) 



First note that the definition above implies that e < L, 
an inequality that will be used below. 

We will define ai and (3i to be some rational approx- 
imations to numbers Inyi and Inxi respectively. Let (jji 
and 9i represent the deviation of these rational approxi- 
mations from the true values, 



ai 



In yi 
In Xj 



(38) 
(39) 



As the rational numbers are dense, these deviations can 
be chosen essentially arbitrarily. But, for our purposes. 
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we are going to choose them as 



7;— < 0i < — for 1 < i < n — 1 
2n n 



< e 



(40) 
(41) 



In other words, the rational approximations ai for aU 
elements excepting the last one are to be chosen such that 
the corresponding deviations (f)i are positive and small, 
but they are also required to be sufhciently far away from 
zero. The last element is an exception. In that case a„ 
has to be chosen as a rational number so that this time 
the sum in (I4ip is made very small. In that case, </>„ 
does not need to be positive. Note that the conditions 
(|40| and ([4T|) provides n separate intervals to choose a; 
from. As rational numbers are dense, all of can be 
chosen as rational numbers. Similarly, we define Pi and 
the corresponding deviations 9i such that 



< 9, < 



2n 



for \ < i < n — 1 



i=l 



(42) 



(43) 



where the deviations for the first n— 1 elements are chosen 
this time to be negative. Similar comments apply in here. 

Below, however, we will need a uniform bound on 
all the deviations. For this purpose, note the following 
bound on (hr, 



n-l 



VnlM < + 2^ 2/^|0^| < + (1 - Vn)- (44) 



1=1 



l-^nl < -+(--1)- 
Vn \yn J n 

< ne^ + (n - 1)- < e 



(45) 
(46) 



where we have used the fact that ?/„ > 1/n for the max- 
imum element of y. Therefore, the following uniform 
bounds can be placed on all deviations 



<e , \9^\<e fori = 1,2, 



(47) 



where the bounds on 9i follow by a similar analysis. For 
most of the following, we will use these uniform bounds. 
The stricter bounds given in ([^0]) and ([1^ will only be 
necessary at the very end. The following bounds on the 
rational approximations will be occasionally used: |Q;i| < 
I Inyil + < L + e < 2L and similarly \/3i\ < 2L. 
Consider the following function 



F(A) = (e^"* - e^^' 



i=l 



(48) 



Our first job is to establish that this function has a root 
near 1, i.e., there is a number Aq, which is very close to 



1 such that -F(Ao) = 0. Once this problem is solved, the 
two new sequences x and y can be defined as 



=.Ao/3i 



(49) 
(50) 



where Zq = YJ^^i e^«"' = X^iLi ■ 1^ that case, both 
X and y are normalized sequences. However, in order to 
reach to the final conclusion, we also need to place bounds 
on the deviation of both Aq and Zq from 1. Therefore, 
the following analysis of bounds is needed. 

First, we must show that F{X) has a root somewhere 
near 1. For this purpose, we look at the value of F{1). 
By using the following inequalities satisfied by the expo- 
nential function, 1 + < < e* < 1 + < + for all |t| < 1, the 
following bounds can be placed on the first term of i^(l). 



i=l i=l 



(51) 

n n 

Y.e'^' <Y,y,{l + ^, + ^1)<l + 2e^ . (53) 



i=l 



Same bounds can also be placed for the second term as 
well, which lead to 



|F(1)I <3e^ 



(54) 



a very small quantity, which indicates that a root is very 
close to 1. 

However, to verify that there is root around 1 and to 
place a bound on the deviation of the root from 1, we 
must make sure that the derivative F'{X) does not rapidly 
go to zero around A = 1. For this purpose, a lower bound 
will be placed on the derivative for |A — 1| < e/L. First, 
note that 



n 
i=l 



E 



yiCti e 



(55) 
(56) 



and the argument of the exponential is small as 

|(A-l)lny, + A0,|<-^i+(l + -^)e<3e . (57) 

Now, using |e* - 1| < |i| < 2\t\ for all \t\ < 1, we can 
find the following lower bound on the expression above 



E' 

i=l 



> -cr(y) - - 2L • 6e (58) 

> -cr(y) - ULe (59) 
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Similar analysis for the second term of F'{X) gives 

n 

'^(3^e^'^' < -a{x)-13Le . (60) 



where the last inequality can be obtained simply by in- 
specting p7p . In other words, we have Xi > Xi for all 
i < n. The conclusion x ^ x then follows. By the same 
method, it can be shown that y is majorized by y as 



Both of these give the following lower bound on the 
derivative F'{X) for |A - 1| < e/L, 



In^ = + (Ao-l)a. -InZo 



(71) 



F'{X) >H~26Le > -H 



(61) 



> — - ( 2 + £2 > , (72) 



e 

2n 



H J 



By using the lower bound given above it is possible to 
see that F{l + e/L) is positive and F{l — e/L) is negative. 
This guarantees the presence of the root in the specified 
interval. But, this interval is too large for our purposes, 
and we need to find a better bound on the place of the 
root. Using F{Xo) = 0, we can get 



in other words yt > yi for all i < n and therefore x -<t y- 
Finally, we have 

D{x; x) = max 161^ - (1 - Ao)/3, - In^ol (73) 

i 

24L' 



< e + ( 2 + — I < eo 



(74) 



F{1) 



F'{X)dX 



\F{1)\ > \Xo-l\- , 
|Ao-l|.™.^ 



(62) 
(63) 
(64) 



In other words, the root is very close to the value 1. 

One final bound, this time a boimd on In Zq will be 
needed. For this, we first note that 



and similarly D{y;y) < eg. Therefore, the inequalities 
([6][8l) are also satisfied by x and y. It is easy to see that 
X and y satisfy the conditions of case A. The number oj 
is given as exp(Ao/A/') where Af is the common denomi- 
nator of the rational numbers and (3i. As a result, the 
conclusion x -<t y follows. Combined with x ^ x and 
y < y, it leads to the desired result x -<t y-O 
Case C. y has zero components. 

Without loss of generality, it is supposed that x and 
y are normalized, they are arranged in non-decreasing 
order and have no common elements. Let y have m zeros, 
i.e., yi = = • • • = ym = and < y„i+i <•••<?/„. 



Zo^Y. e""e(^°-i)"- < 51 ^"j 6+''^°"''"^ , (65) Let be a sequence defined as follows 



1=1 \i=l / 

and a similar analysis for the lower bound gives 



e for « = 1, 2, . . . , TO 

(1 — me)yi for i — m 



1,. 



. ,71 



(75) 
(76) 



llnZol < 
Finally, ^ and ^ gives 



2|Ao-l|L 



< 2e^ 



(66) 



(67) 



where we have used the fact that i — 1 > Ini > {t — l)/t. 
As a result, we get 



InZnl < 2 



12L\ 



(68) 



Now, it is possible to show that the sequences x and y 
satisfy all the required properties to complete the proof. 
First, we will show that x is majorized by x. For this 
reason, we will look at the ratio Xi/xi ior i — 1,2, . . . ,n— 
1, i.e., for all elements except the last one. Here, we will 
make use of the upper bounds given in (|42p as 



where e is a non-negative parameter. We will only be 
interested in the values of e in the range e < (y.^+i + 
m)~^ , where z*^ is arranged in increasing order. It is easy 
to see that all such sequences are related to each other by 
the majorization relation, i.e., if > es then z"^^ -< z*^^. 
As z^ = y, we have z'^ ^ y for all values of e in the range 
considered. 

Our job is to show that if e is sufficiently small, then 
X and satisfy the inequalities (IB][5]). This is a straight- 
forward but laborious procedure which is detailed below. 
For this purpose, different intervals of v values will be 
considered separately and for each interval, the existence 
of a separate upper bound for e will be provided. 

(a) For v < 0; The quantity ei = yn{xi/yn)"'^"^ is a 
possible upper bound for this range. Let e < ei. For the 
special case v = Q,we have 



Mx)_ 

Ao(z') 



> — - 2 
- 2n 



+ (l-Ao)/?, +lnZo 
24i 



H 



> 



(69) 
(70) 



> 



e"(l - me)"- 
2/ri ^ e 



nn 
i=m+l Vi 



(77) 
(78) 



For all negative values of i' we make use of Bernoulli's 
inequality, which states that a'' — 1 > r(a — 1) for any 
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r > 1 and any positive number a, to reach 

mle-^ - > - • (79) 
This then leads to 



i=l 



^(z^" - me'^ + il^mer ^ y." (80) 

i— m+l 

n 

> me" + {n — m)y'^ > nx^ > xj' (81) 



i=l 



As a result, we conclude that A^{x) > A^{z'^) for all 
u <0 whenever e < ei. 

(b) For Q <u <l/2: The function 



Ju 



(82) 



is strictly positive in the interval (0,1/2] and moreover 
it has a strictly positive limit at = 0. Therefore, £2 — 
min^g[o.i/2] Jv is a positive number. If e < £2, we have 

n n n 

^xr>m£^+ ^ 2/r>E«r ' (83) 

i— 1 i— m+1 i—1 

which leads to A^{x) > A^{z'^) in this interval, 
(c) For 2 <v: Let if be defined as 



K = max 



Au{x) 



(84) 



which is a positive number such that K < 1. Note that, 
as X and y have no common elements, the ratio above 
dX V ^ +00 gives x\/y\ which is smaller than 1. Let 
£3 = (1 — K)/m. Then, for any £ < £3 and for all > 2 
we have 



Y,{zt)^ > [l-meY yr 

n n 



i=l 



(85) 
(86) 



i—m+l i—1 

This shows the desired inequahty, A„{z'^) > Ai,{x). 
(d) For l/2<iy< 2; Let 

i?.=^in444 . 

u-1 A^x) 



(87) 



The inequalities ([B][5]) imply that R,^ is a strictly positive 
continuous function in the interval considered. Therefore, 
the minimum M = min^g [1/2.2] is a positive number. 
Let 

Since all sequences z'^ are related into each other by the 
majorization relation, for any eA > £b we have R^{eA) < 



Rv{^b) for all V. In other words, as £ decreases, the 
function R^{e) monotonically increases. Finally, we note 
that Ru{() converges pointwise to R^ as £ goes to zero. 
At this point, we invoke Dini's theorem, which states 
that a sequence of monotonically increasing, continuous 
and pointwise convergent functions on a compact space 
are uniformly convergent. Therefore, there is a positive 
number £4 such that whenever £ < £4, we have Rv{e) > 
M/2. 

For such values of £, the inequalities ^ and ([7]) are 
satisfied for all v £ [1/2,2]. Moreover, the inequality ([8]) 
is also satisfied, as Ri{e) — a-{x) — o-{z'^) > M/2 > 0. 

As a result, if £ < min(£i, €2, €3, €4), then the sequences 
X and z"^ satisfies all the inequalities ([B][5]). The proof of 
case B enables us to conclude that x -<t ■ Finally, by 
z^ ^ y we reach to the desired result x V- 



IV. DISCUSSION AND CONCLUSION 

A set of necessary and sufficient conditions are given 
for the trumping relation. The conditions involve a con- 
tinuous variable, but they are easy to verify for concrete 
examples. 

Conditions can be easily adopted to the sequences 
in the closure of T{y), where 



T{y) = {x -.x y} 



(89) 



is the set of sequences trumped hy y. In that case, if 
X e T{y) then the conditions ([MH]) must be satisfied but 
with strict inequalities replaced with non-strict ones. 

If X G T{y), but X is not trumped by y, it means 
that no catalyst can achieve the conversion of x into y 
with probability 1, but it is possible to find a sequence of 
catalysts (with growingly large Schmidt numbers) such 
that the conversion probability is made to approach 1. 
Interestingly, this property is also shared by states that 
are far from the boundary of T{y). Consider the example, 



2 3 4 
9' 9'9 
1 2 2 

5' 5' 5 



(90) 
(91) 



As x\ > yl, X is not in the closure of T{y). However, it 
can be verified that A^{x) > A^{y) for all u < 1. This 
then implies that, any given probability less than 1 can be 
achieved by a suitable catalyst in the conversion of x into 
2/.[l3] However, the elements ofT{y) satisfy an additional 
property, i.e., they can be catalytically converted with 
unit probability to another state only slightly different 
from y. It is puzzling to see that this property is not 
shared by the pair x, y given in the example above. 

Once it is understood that catalysis is possible, the 
problem of finding a suitable catalyst can in principle be 
solved by going backwards along the proofs. Although 
possible solutions of the problem posed in the Lemma 
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in Section |TT] can be found by the well-established pro- 
cedures of linear programming, carrying out the whole 
procedure for realistic cases might be forbidding, as the 
degree of the polynomial 7(5) and of the sought for poly- 
nomial a(s) might be very large. However, the method 
used in the proof of the Lemma can used to place an up- 
per bound on the degree of a{s) (but not on the Schmidt 
number). This also suggests a conjecture that the com- 
plex roots, I/, of the equation A„{x) = A,^{y), and their 
closeness to the real line could be used for estimating the 
minimum amount of resources the catalysts should have. 



APPENDIX A: PROOF OF THEOREM 2 

The most troublesome part of the proof of Theorem 2 
is the neighborhood of j/ = 1. This part can be handled 
with the following theorem. 

Theorem 3. For any positive sequence x and any 
given S > 0, there is a positive number e such that, for 
any x with x,; — '^Xi, and D{x; x) < e we have 

e-^l"--!! < 44^ < e'l-^-^l Vi/e [1/2,2] . (Al) 
Au(x) 

Proof of Theorem 3. Without loss of generality, it is 
supposed that x — x"^ . Let S,^{x) — be the i^th 

power sum and Ki, be defined as 



K„ = In 



S„{x) 
S„{x) 



(A2) 



Note that Ki = 0. We place the following bound on the 
absolute value of v derivative of K^, 



dv 



1 



< e 



< e 



Sy{x)Su{x) 
1 



Su{x)S^{x) 
1 



Su{x)S^{x) 



5:<xpn|i (A3) 

E^r-P-S (A4) 

Y,{xtx';-x';x-)\n- 



X-i 



Since, for i > j we have Xi > Xj, all of the logarithmic 
terms are non-negative in the expression above. As a 
result, for any positive v, 



dK,. 



dv 



< e- 



< e 



< e 



S„{x)Si,{x) 



^<xpn^ (A5) 



2Su{x)Su{x 



EV V 
Xj Xj 



1 '^i 

in — - 



e "'^ Siy{x) Xn 
■ in — 



S^{x) xi 



(A6) 
(A7) 



Finally, we can apply S,y{x) > e '^'^S„{x) to the last line 
which gives 



dK^ 



dv 



Jlvt 



< e 



< e 



1 1 '^n 

In — 

2 xi 

2 xi 



(A8) 
(A9) 



where the last inequality is valid for all < < 2. Note 
that the right-hand side of the last expression has zero 
limit as e — > 0. This enables us to choose the value of 
e so small that the right hand side is less than 6/2. In 
other words, \dK^/dv\ < 6/2. 
Next, we express K^, as 



" dK„ 



dv 



-dv 



The inequality above then implies that 



1 



\K,.\<-6\v- 



(AlO) 



(All) 



Finally, considering only the values of v in the interval 
[1/2, 2], we have 



In 



A^{x) 



A^{x) 



\M<^Jt^<s\v- 



(A12) 



which is the desired result. □ 

Proof of Theorem 2: Without loss of generality sup- 
pose that X and y are normalized, i.e., '^Xi = — 1. 
Since the minimum and maximum values of the these 
sequences are different by the assumptions of the theo- 
rem, the strict inequalities are valid at the infinities, i.e., 
A-ooix) > A-oc{y) and Aooix) < Aoo{y). Let d, be 
defined as 



G„ = In 



A^{x) 



(A13) 



By the inequalities (IB][5]), we have < for all < 1 
and > for all > 1. At infinities approaches to 
non-zero limits. Moreover, the derivative of a,t v ~ 1 
is 



G[ = a{x) - a{y) > 



(AM) 



Therefore, both of the following quantities are strictly 
positive, 



B = min \Gi, 

i/e[-oo,l/2]u[2,oo] 



M 



G, 

mm 

^£[1/2, 2] V — 1 



(A15) 
(A16) 



By Theorem 3, there are numbers ei and £2 such that 
D{x;x) < ei implies that \lnA^{x)/A^{x)\ < M\v-l\/3 
and D{y]y) < €2 implies that | ln>l^(y)/yly(j/)| < M\v — 
l|/3. We choose e = min(ei, £2, i?/3). 



Let X and y be arbitrary sequences such that ^ Xi 
Y,m = l, D{x\x) < e and D{y;y) < e. Let 



Gi/ — In — — ^^-J- — Gi, 
Ay{x) 



Ay{y) A^{x) 



Our purpose is to show that satisfies the desired prop- 
erties, i.e., it is negative for z/ < 1, positive for v > I and 
has a simple zero at ly ~ 1. Note that D{x; x) < e implies 
that I \nA^{x)/A„{x)\ < e for all v. 

We consider the following ranges of v values separately, 

(a) For iy<l/2, we have G^ < -B + 2e < -B/3 < 0. 

(b) For iy>2,we have G^ > B -2e> B/3 > 0. 

(c) For v e [1/2, 2] we have 
G„ , . 2M 



1 



> M 



> 



As a result, Gj, satisfies the desired properties in 
this interval as well. 



Finally, for the inequality ([5]), we note that for e 
[1/2, 1), we have 



G. > ^ 
- 1 - 3 



(A19) 



Taking v I limit gives the derivative of G^ which is 



G[=a{x)-a{y)>^ 



(A20) 



(A18) This completes the proof of Theorem 2. 
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